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^~^' Abstract. We introduce our package +Ideals for Magma, designed to per- 

C^ ' form the basic tasks related to ideals in number fields without pre-computing 

integral bases. It is based on Montes algorithm and a number of local tech- 
niques that we have developed in a series of papers in the last years. 

in 

(N 

C~' ' Introduction 

^^ , A commonplace in Number Theory problems is the need of effective computa- 

rS^ • tions. The last decade has seen an exponential growth of mathematical software 

j^ I to satisfy this necessity. Striking ideas have expanded our computational limits far 

beyond our present requirements in many fields, specially those related to crypto- 
graphic applications. However, there remain some parts of Computational Number 
Theory which have not evolved that much. Algebraic Number Theory appears to 
be one of these. 

While the available algorithms for working over number fields are bright and 
O^ I efficient, the computational requirements of real life problems in Number Theory 

in ' frequently exceed their capability. As soon as we have to work in a number field 

\l I defined by a polynomial with large discriminant, the determination of its ring of 

l/^ ■ integers, which is the cornerstone of all the classic algorithms, becomes unfeasible, 

^^ I thus making impossible further computations. 

^^ ■ We have developed a number of techniques, based on higher order Newton 

polygons, to deal with the prime ideals in number fields ( |GMN08j . |GMN09a] . 
[GMNOQb] ■ [GMNIO) ). The application of these ideas leads to extremely fast and 
is^ I efficient algorithms for the basic tasks in Computational Number Theory. These 

'j_j ■ new algorithms can work over number fields of degree 1000 or with a 300-digit 

C^ I discriminant with a home computer in a few seconds. A new philosophy is beyond 

this project: the centre of attention in Computational Algebraic Number Theory 
may be not the ring of integers, but the prime ideals. 

We have implemented these algorithms in the +Ideals package for Magma. 
Besides factorization of ideals, this package gives facilities to compute valuations 
at prime ideals, reduction modulo prime ideals and Chinese Remainder Theorem 
problems. The algorithms have been presented in I GMNIOJ . The package can 
also compute integral bases of number fields, following the algorithm introduced 
in |GMN09a] . The goal of this paper is to introduce the package, presenting some 
examples of use and describing its main functions. 
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The paper is structured in three sections: after giving a general overview of the 
package in the first section, we illustrate its use and capabilities by means of some 
concrete examples in Section 2. Section 3 gives a more exhaustive description of 
routines of the package. 

A primitive version of the package was published about two years ago, under 
the name Newton. It included essentially the decomposition of primes in number 
fields and the computation of integral bases. The functions in the Newton package 
remain in the +Ideals package, to ensure compatibility for all users. 

The +Ideals package can be obtained on demand from the authors, or down- 
loaded from its web page ^ http: //www-ma4.upc . edu/~guardia/+Ideals .html) . 
The package requires version v2.15 of Magma or newer. 

1. General overview of the package 

The +Ideals package is concerned with the arithmetic of fractional ideals in 
number fields. The keystone is Montes algorithm, which provides local information 
for the decomposition of rational primes in number fields. The great advantage of 
the algorithm is to avoid the use of integral bases. This shortcut makes possible to 
perform fast computations in number fields of huge degree and index. 

The underlying philosophy of the package is that prime ideals are the basic tools 
for many computations in number fields, and hence all the functions in the package 
perform tasks related to them. The main task, of course, is factorization of frac- 
tional ideals, but computation of valuations and reduction maps are also important. 
As soon as a prime ideal is detected along a computation, the information which 
can be useful for other computations is stored in some new attributes for number 
fields, so that it can be easily retrieved if needed. 

The package introduces two new data types, PrimeldealRecord and Ideal- 
Record, designed to represent fractional ideals in number fields, and some attributes 
for number fields to store information about these ideals. Most of the functions 
included in the package are extended versions of the usual Magma routines for ideal 
arithmetic, designed to deal with these new data types, so that the experienced user 
of Magma should find no difficulty in the use of the package. 

Let us finish this overview with a self-explanatory small piece of code illustrating 
how to perform the most common tasks in ideal arithmetic with our package: 

> Zx<x> :=PolynomialRing(Integers() ) ; 

> K:=NuinberField(x'100-x~75+x~50+2~500) ; 

> I:=ideal(K,2); 

> time Factorization(^I) ; 
Time: 1.300 

>1~ Factor izationString; 
P(2,1)-2*P(2,2)-2*P(2,3)-2*P(2,4)*P(2,5)*P(2,6)*P(2,7)*P(2,8) 

> K'Locallndex[2] ; 
12250 
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2. Some hard examples 

We illustrate the use of the package by means of some exam.ples picked from, 
different topics in Computational Number Theory. The computations have been 
made with Magma v2. 15-11 in a Linux server, with two Intel Quad Core processors, 
running at 3.0 Ghz, with 32Gb of RAM memory. All the running times mentioned 
in these examples are expressed in seconds. 

2.1. A real life problem. A very common problem in Arithmetic Geometry is the 
determination of good models for curves over number fields. Usually, a curve with 
potentially good reduction acquires good reduction over a number field generated by 
certain torsion points in the jacobian of the curve. The number fields encountered in 
this process are usually huge, and it is a difficult computational problem to find the 
places over the rational primes we are interested in. We show in this example how 
to use our package to solve such a problem. We will determine the reduction type 
of a rational elliptic curve over a huge number field, applying Tate's algorithm. The 
decomposition of primes in the number field is the minimal information required, 
but no integral basis is necessary to run the algorithm, so that the problem fits 
very well the philosophy of our package. 
Consider the elliptic curve 

E -.y^ ^x^ + 3x^ + 3a;, 

and let K = Q{0) be the number field determined by a root 6* G Q of the 17-th 
division polynomial. It is a number field of degree 144. We can determine it with 
the following Magma code: 

> E:=EllipticCurve( [0,3,0,3,0] ) ; 

> f :=DivisionPolynomial(E, 17) ; 

> K<w>:=NumberField(f ) ; 

The discriminant of the polynomial / defining K is huge, but easy to factor, due 
to its geometric origin: 

> d: =Discriminajit (f ) ; 

> Factorization(d) ; 

[ <2, 13728>, <3, 10296>, <17, 143> ] 

We cannot expect to find the ring of integers Zk of such a big number field with 
the standard Magma routines: the function MaximalOrder exhausts the com- 
puter's memory after some hours of work. But we do can compute in this ring! We 
must redefine K through a monic polynomial with integer coefficients; since the 
leading coefficient of / is 17, we simply type: 

> Qx<x> :=PolynoinialRing(Rationals() ) ; 

> g:=17'143*Evaluate(f ,x/17) ; 

> K<w> :=NumberField(PolynomialRing(Integers() ) !g); 

Note that a root of g is w — 179. Let us now attach our packagqj: 

> Attach("dir/+ideals.m") ; 



Replace dir with the name of the directory where you saved the package. 
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Let US find the decomposition of 2 in "Lk- We first define tlie ideal 2'Lk using tlie 
new function ideal, and then we factor it into prime ideals: 

> i:=ideal(K,2) ; 

> time Factorization(^i) ; 
Time: 0.010 

> i~FactorizationString; 
P(2,l)~3*P(2,2)-3*P(2,3)-3*P(2,4)-3*P(2,5)-3*P(2,6)-3 

The result means that 22^ is the cube of the product of six different prime ideals: 
2Z7< = (pi • ■ -PqY- The use of ^ allows the factorization routine to append infor- 
mation to the record corresponding to the given ideal. A shorter way to find the 
decomposition of a prime in Z^ is: 

> time Factorization(ideal(K,3) ) ; 

[ [ 3, 1, 2 ], [ 3, 2, 2 ], [ 3, 3, 2 ], [ 3, 4, 2 ] , [ 3, 5, 

2 ], ] 

Time: 0.010 

which indicates that there are 5 different prime ideals in "Lk over 3, all of them 
with ramification index 2 (but now we can not access the assigned attributes of the 
ideal, since they are not saved) We can access these prime ideals through the new 
attribute Primeldeals established for number fields. The command: 

> K~ Primeldeals [2] ; 

will generate a long output with all the data already computed for these ideals. For 
instance: 

> K~ Primeldeals [2 , 1] ' f ; 
8 

gives the residual degree of the first of the prime ideals dividing 2. The residual 
degrees of the prime ideals above 3 are: 

> [P~f:P in K 'Primeldeals [3] ] ; 
[ 8, 16, 16, 16, 16 ] 

We already have generators for the ideals at hand; typing: 

> alpha :=K~ Primeldeals [3, 1] ~ Generator ; alpha; 
1/3433683820292512484657849089281* (w-129+2451*w- 128+ . . . 

we obtain an element a 6 "Lk such that p ~ (3, a) is the prime over 3 with residual 
degree 8. 

Along the factorization routine of the ideal "pliKi the program attaches to each 
prime ideal overp a local ij/pe, which carries on some canonical invariants ([GMNOSb]). 
These types are stored as an attribute of the prime ideals (K~ Primeldeals [p,k] ~Type) 
and contain essential information. As a by-product of the computation of these 
types we obtain the p-index of the polynomial g. It is saved in an associative ar- 
ray of Magma called Locallndex, with one entry for each prime number already 
factored in Z^f : 
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>K~LocalIndex[2] ; 
6816 

Using this array, we can easily find the discriminant of K\ we first factor completely 
the discriminant d oi g and then get the global index: 

> time Factorization(id.eal(K, 17) ) ; 

Time: 0.380 

> ind : = (2' (K~ Locallndex [2] ) *3~ (K~ Locallndex [3] ) * 
17" (K~ Locallndex [17] ) ) ; 

> disc:=d div ind~2; 

> Factorization(disc) ; 

> [ <2, 96>, <3, 72>, <17, 143>] 



We can check that the point P = (w/17, v^zi^Vr73T3zi^Vl72TWl7) G E{K) 
is not defined over K. The computation of the square root in the y-coordinate may 
be quite expensive, so that we simply reduce its second coordinate modulo a prime 
ideal: 

> Factorization(ideal(K,5)) ; 

> yred: = (w'3/17'3+3*w'2/17'2+3*w/17) mod K~PrimeIdeals [5, 1] ; 

> IsSquare(yred) ; 
false 

Let us apply now Tate's algorithm to determine the reduction type of E over the 
prime p — (3, a) defined above. We have ai = as = qq = 0, 02 = 04 = 3, 62 = 12, 
64 = 6, fee = 0, ^8 = ^9. It's a simple matter to check the cumulative hypothesis of 
Tate's algorithm: 

> P:=K~PrimeIdeals[3,l] ; 

> dE: =Discriminant (E) ; 

> Evaluation (dE,P) ; 
6 

> [Evaluation (ai, P) : ai in alnvariants(E)] ; 
[ Infinity, 2, Infinity, 2, Infinity ] 

> [Evaluation (bi, P) : bi in blnvariants(E)] ; 
[2, 2, Infinity, 4 ] 

Thus, we have 

t'p(A) >0 

Vpias) > 0,Wp(a4) > 0,Wp(a6) > 0, 

Vpib2) >0, 

Vp{ae) > 2, 

Vpih) > 3, 

Vpibe) > 3, 

Vp{ai) > 0,Wp(a2) > 0,Wp(a3) > 2,Wp(a4) > 2,Wp(a6) > 3. 

We must now consider the polynomial 

P(T) = T^ + a2,iT^ + a4.2T + a6,3, 

where ai^m = Oj/n™, and 11 is a local uniformizer for p. We get one by typing: 
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> PI :=P~ Generator; 

Since 02,1 = 05,3 = O(niodp) and Wp(a4.2) = 0, the polynomial P(T)(modp) has 
three different roots in an algebraic closure of the residue field k{p) of p, so we 
already know that the Kodaira symbol of i? is Jq . The local conductor is p"t)(^)"'' = 
p^. To determine the local index Cp = [E{Kp)/Eo{Kp)], we need the number of 
roots of P(T) in fc(p). We may proceed as foUowiO : 

> time m:=3/PI'2 mod P; 
Time: 13.570 

> F:=ResidueField(P) ; 

> FT<T>:=PolynomialRing(F) ; 

> Factorization(T~3+m*T) ; 
[ <T, 1>, 

<T~2 + 2*zl'7 + 2*zl'6 + 2*zl~4 + 2*zl~3 + zl + 2, 1> ] 

Hence, the local index is Cp — 2. It is easy to check that for the remaining prime 
ideals in Z^ above 3, the algorithm proceeds in the same way (all the ramification 
indices are 2) up to this point, where we find that the corresponding polynomial 
P{T) has three different roots in the residue field, and hence the local indices at 
these primes are equal to 4. 

Finally, our package has its own routine IntegralBasis to compute a Z-basis 
ofZic: 

> time basis :=IntegralBasis(K) ; 
Time: 18.620 

It has to be mentioned, however, that this routine is based on a conjecture that 
we have not been able to prove yet f |GMN09"a] l. In any case, the routine always 
returns a correct result (it performs a check of its output) and is prepared to warn 
the user in the improbable case that the algorithm fails (after millions of tests, we 
have not encountered any counterexamples). 

2.2. A tall number field. In this second example we will force the computational 
capabilities of the package to show its power. The concrete computations are not 
significant on its own; they have been chosen to illustrate the different routines of 
the package. 

Consider the number field K — Q{6) given by a root G Q of the polynomial 
f{x) = x^"™ + 2-'^°a;5° + 2^°. The factorization of the discriminant of / is 

Disc(/) = 2539^°3^°52™0i27503i350743504ggg22952750p50^ 

with p = 337572698551220494882323528404563236947916489629537. The large de- 
gree of / makes impossible to work in this number field using the standard func- 
tions of Magma, even after factorizing the discriminant, since the computation of 
the integral basis is compulsory for these functions. But our functions avoid this 



The computational representation of the residue field of a prime ideal is computed as a 
successive extension of finite fields: 

Fo CFi C ■■■ CF,.+i =Fp =¥olzo,zi,...,Zr], 

where Fq is the prime field and F^ := Fi_i(2i_i), for f < i < r + 1- The parameter r is the order 
of the type attached to p. 



ARITHMETIC IN BIG NUMBER FIELDS: 



THE •+IDEALS' PACKAGE 



computation, so that we can work with ideals in K. Let us begin by determining 
the discriminant of K: 

> Attach (" /dir/+Ideals. m" ) ; 

> Z:=Iiitegers() ; 

> Zx<x> :=PolyiiomialRing(Integers() ) ; 

> f := x'lOOO + 2'50* x'50 + 2~60; 

> d:=Discriminant (f ) ; 

> Factorization(d) ; 

[ <2, 53940>, <3, 50>, <5, 2000>, <127, 50>, 
<313, 50>, <743, 50>, <4886229527, 50>, 
<337572698551220494882323528404563236947916489629537, 50> ] 

> badprimes :=PrimeDivisors(d) ; 

> K<w>:=NumberField(f) ; 

> for p in badprimes do time Factorization(ideal(K,p) ) ; end 
for; 



> [K~LocalIndex[p] : p in K~ Factor izedPr imes] ; 
[ 26235, 0, 20, 0, 0, 0, 0, ] 

The last command uses the new attribute FactorizedPrimes for number fields: 
it is a list containing the rational prime numbers that have already been decomposed 
in Zk already. The running time of the factorizations above and the local indexes 
are gathered in the following table: 



Ideal 


Index 


Time 


2Zk 


26235 


0.36s 


3Zk 





0.61s 


5Zk 


20 


0.63s 


127Zk 





1.29s 


313Zk 





3.69s 


743Zx 





6.47s 


4886229527Zk 





6.96s 


PZ^K 





60s 



Thus, we need less than 90 seconds to see that the discriminant of K is 
Disc(ii') === 2i'*™3^°5i9«°127^°313''^"743^°4886229527''^V°- 

While we cannot expect to factor the ideals / = (ui^ + 50)Zx, J = (w + 10)Zk 
in a reasonable time, because this implies the factorization of their norm, we can 
factor their sum in a moment: 

> I:=ideal(K,w'3+50) ; 

> J:=ideal(K,w+10) ; 

> N:=I+J; 

> time Factorization(~N) ; 
Time : . 030 

> N~FactorizationString; 
P(2,1)-2*P(2,2)-2*P(2,3)~2*P(2,4)-2*P(2,5)'2 
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This output means that ah the primes P(2,l),..., P(2,5) appearing in the 
hst K'Primeldeals [2] divide I + J with exponent 2. Note that, besides defining 
them with the new funetion ideal, we manipulate the ideals with the standard 
Magma notations. We can also multiply and divide fractional ideals in K. 

> R:=ideal(K,2); 

> Factorization('-^R) ; 

> time S:=R/N; 
Time: 0.280 

> S~FactorizationString; 
P(2,l)~8*P(2,2)-8*P(2,3)-36*P(2,4)-36*P(2,5)-36 

Generators of integral ideals are efficiently computed, but the computation of 
generators of non-integral ideals could be quite time-costing. In the above example, 
it would take much time to determine generators for N^'^, but: 

> time gens :=TwoElement (S) ; 
Time: 0.000 

Our last computation in this field in this number field will be a p-adic square 
root. We take the first two prime ideals pi,p2 in the list K~PrimeIdeals [5] , and 
check whether w/2 is a square in their residue fields; 

> Pl:=K~PrimeIdeals[5,l] ; 

> P2:=K~PrimeIdeals[5,2] ; 

> IsSquare(w/2 mod PI); 
false 

> t,xO:=IsSquare(w/2 mod P2);t;x0; 
true 

2*zl+4 

After Hensel's lemma, we now that w/2 has a square root in Kp^, but not in 
Kp-^. We can compute this root with Newton's method, but we will need to lift 
elements from the residue field of pi to K. This task is performed by the routine 
Lift. For instance, to lift xO we make: 

> X:=Lift(xO,P2); 

The following code computes \fwj2 G Kp^(mo'i)(i^y. 

for n in [1. .10] do 

lambda :=Lift((2*X mod P2) ' (-1) ,P2) ; 

X : =X- (X~2-w/2) *lambda; 
end for; 
print PValuation(X'2-w/2,P2) ; 

The whole computation takes about 60 seconds. 



2.3. A fat number field coming from the modular w^orld. The degree of a 
number field has a strong influence in the cost of the computations on it, but it is 
not the only handicap, nor the principal one. The discriminant turns out to be the 
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highest obstacle in many cases, since it may be a very large integer, impossible to 
factor. The following example illustrates this situation. 

Let / G ro(l) be the weight 76 modular form given by the Magma command 

> f := Newf orms(CuspidalSubspace(ModularForms(l,76))) [1, 1] ; 



It is defined over the number field K = Q{9), where 6* £ Q is a root of the 
polynomial 

F{x) = x^+ 57080822040x5 - 198007918566571424544768a;'' 
-11405115067164354385292006554337280a:3 
+9757628454131691442128845013041495838774263808x2 
+2900139955623795004984359750037160248001145937615808102402: 
-92217203874207784163935379997152082331434364841943058919508374716416. 

We can access this number field in Magma with the command: 

> K := ParentCCoefficientCf ,1)) ; 
The discriminant of F{x) is 

disc(F) = 22'^''3^25^^7^^112i3217^19243259-1932293-3919872475942721376796815212m, 

where M is a composite integer of 135 decimal figures which we have not been able 
to factorize. J. Rasmussen asked us (fRasm ussenlO| ) for a test to check certain 
divisibility conditions on the ring of integers of the number field K — Q{0), related 
to his work on congruences satisfied by the coefficients of certain modular forms. 

Since we cannot factor completely the discriminant, the Magma function MaximalOrder 
is useless. But we can do many computations in K with our package. Let us per- 
form a Chinese Remainder computation. We first define the number field and find 
the decomposition of 3 in the ring of integers: 

> Attach("dir/+ideals.m") ; 

> I:=ideal(K,3); 

> Factorization(^I) ; 

> I~FactorizationString; 

P (3 , 1) *P (3 , 2) *P (3 , 3) *P (3 , 4) *P (3 , 5) 

The factorization routine automatically finds generators for the prime ideals it 
encounters. They are stored in the attributes IntegralGenerator and Generator 
of the IdealRecord record subjacent to the ideals, but they can be accessed with 
the TwoElement function: 

> TwoElement (K~ Prime Ideals [3,1]); 
[ 3, 

1/177147* (2*w'5 + 1815*w"4 + 586980*w~3 + 732159*w"2 + 
658287*w + 1535274) ] 

Consider the following Chinese Remainder problem: find an element a (z K 
satisfying 

a = w(modpi), a = w"^ {mod p"^) , a = w^(modp3), 

a = i(;'*(modp|), a = l(modp5). 
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where pi,...,p5 are the prime ideals above 3 in Zx, ordered according to the 
output of our factorization function. We can solve the problem using our extended 
version of the standard CRT function of Magma. This extended version takes a list 
of numbers in a number field and a list of ideals in its ring of integers, and solves 
the corresponding congruences. The following piece of code finds a solution to our 
problem: 

>ids : = [K~ Primeldeals [3 , 1] , K~ Primeldeals [3 , 2] ~2 , 

K~ Primeldeals [3 , 3] "3 , K~ Primeldeals [3 , 4] ~4 , K~ Primeldeals [3,5]]; 

>time r:=CRT([w,w"2,w"3,w~4,l] ,ids) ;r; 

Time: 0.010 

1/19683* (1132037*w-5 + 1583958*w-4 + 498663*w"3 + 923157*w'2 + 

649539*w + 354294) 

We can check the result using the PValuation routine: 

>PValuation(r-w,K~PrimeIdeals[3,l]) ; 

7 

>PValuat ion (r-w~2 , K~ Primeldeals [3,2]); 

4 

>PValuation(r-w~3,K~Primeldeals [3,3] ) ; 

4 

>PValuat ion (r-w~4 , K~ Primeldeals [3,4]); 

4 

>PValuation(r-l ,K~ Primeldeals [3,5]); 

1 

A second partial test of the result can be done reducing the element r modulo 
pi and p2: 

>(r-w) mod K 'Primeldeals [3, 1] ; 



>(r-l) mod K'Primeldeals[3,5] ; 



2.4. A two-parametric family of examples. We now present a two-parametric 
family of polynomials, which constitutes a good bench mark for our package, since 
the parameters involved give the chance to grade the difficulty of the computations 
at our convenience. 

Let 0„ G Q be a root of the polynomial (?„ = a;" + a; — 1, and let Kn = Q(^n)- 
We assume that gn is irreducible. For every r G Z, the element 6n,r '■= '''On is also 
a primitive element of K^ with minimal polynomial /^ = x" + r^^^x — r". The 
discriminants of these polynomials are: 

Disc(5„) = (_l)(»-l)(«-2)/2(„« ^ (^ _ 1)"-^)), 

Disc(/;) = r"("-i) Disc(g„). 

Suppose we are given /^ as defining polynomial for K^. Those primes dividing r 
but not A((7„) will not ramify in Kn and we can make the local index arbitrarily 
large for them. Thus, an (im)proper election of r can make computations in the 
field Kr, as difficult as we want. 
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We will use the following function to make some tests: 

> function test(n,r) 

function> K:=NumberField(x~n+r~ (n-1) *x+r~n: Check:=f alse) ; 

f unction> a:=Factorization(ideal(K,r) ) ; 

function> return [K~LocalIndex [p] : p in PrimeDivisors(r)] ; 

function> end function; 

The function takes the degree n and the multiplier r, and determines the decom- 
position of the primes dividing r in the field Kn, defined by means of the polynomial 
/^, and it returns the local index of the discriminant of /^ at them. The following 
table collects the running times of the test for some values of n and r: 



r\n 


100 


201 


300 


400 


501 


250 


0.02 


0.05 


0.33 


1.04 


0.69 


2iUU 


0.05 


0.1 


0.85 


2.68 


1.77 


2^00 


0.07 


0.28 


2.23 


7.4 


4.73 


2400 


0.17 


0.76 


6.15 


20 


13.71 


650 


0.13 


0.77 


1.7 


2.9 


9.77 


glOo 


0.31 


2.04 


4.57 


7.75 


28.09 


g^OO 


0.5 


5.71 


12.78 


21.75 


78.38 


30^" 


0.35 


1.49 


5.02 


6.980 


20.910 


30100 


0.82 


3.91 


13.7 


18.34 


58.81 


302O0 


2.24 


10.95 


38.49 


49.66 


630.81 


210^" 


0.64 


2.31 


9.63 


10.06 


35.25 


210^"" 


1.48 


6.03 


26.62 


25.98 


99.41 


2310"'" 


2.63 


12.25 


31.72 


44.33 


161.55 



For r — mJ' and fixed n, the local indices at the primes dividing ra are all equal to 
kn{n — l)/2. For instance, the largest local index in the table was 50100000, in the 
case k = 400, n = 501. 

We can make some comparisons with the standard Magma routines for number 
fields. For instance, we have compared the running time of Magma's command 
MaximalOrder and our routine IntegralBasis, with the following results|j: 



These running times do not include the time spent factoring the discriminant of the polynomial 
/^; which is performed before the computation of the maximal order. 
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n 


r 


MaximalDrder 


IntegralBasis 


5 


2t>u 


0.06 


0.01 


5 


2 100 


0.25 


0.02 


5 


2^00 


1.77 


0.12 


5 


2400 


19.58 


0.84 


5 


gbU 


0.55 


0.04 


5 


glUU 


2.81 


0.23 


5 


6^00 


23.76 


0.01 


5 


64UU 


270 


0.02 


10 


250 


1.6 


0.18 


10 


2 100 


14.37 


1.11 


10 


2^00 


185 


0.1 


20 


250 


21.37 


0.06 


20 


2 100 


185 


0.1 


40 


210^ 


168 


0.71 


40 


210^^ 


1066 


3.28 



One may wonder whether Magma's factorization routine is faster than ours if 
the maximal order has aheady been computed. The following table compares the 
performance of both routines for n = 40 and muhipliers ri ~ 210^ and r2 = 210^: 



r = 210"^ 




Magma 


+Ideals 


2Zk 


0.35 


0.02 


3Zk 


1.77 


0.00 


5Zk 


0.91 


0.01 


7I^K 


1.11 


0.01 



r = 210^ 




Magma 


+Ideals 


2Zk 


0.36 


0.01 


3Za- 


1.78 


0.01 


5Za' 


0.92 


0.00 


7Zk 


1.16 


0.01 



The local index at these primes is 2340 for r = 210^ and 3900 for r = 210^ 



2.5. The other way round: how to construct hard examples. This last 
example illustrates the reverse use of the package to construct examples of number 
fields with prescribed decomposition of a given set of primes, determining the fiux 
of the Monies algorithm for them. This section may be skipped for the reader not 
specially interested in the theory behind the package. 

According to the philosophy introduced in |GMN10] , a prime ideal p in a number 
field can be canonically represented as a list 



b; 



^pj: 



where the (pk are certain monic polynomials with integral coefficients, determined 
by Montes algorithm, when applied to the polynomial defining the number field and 
the rational prime p. These polynomials are the main data in the type corresponding 
to the ideal p; the parameter r is the order of this type. 

Let us explain how to build a type t with the package. We fix the prime p = 2 
and the irreducible polynomial (pi = y + 1 £ F2 [y] , and initialize our type 
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> Attach (" /dir/+Ideals. m" ) ; 

> ZX<x> :=PolynomialRing(Integers() ) ; 

> FY<y>:= PolynomialRing(GF(2) ) ; 
t ,Y,z:=IiiitializeType(2,y+l) ; 

The variable t will contain our type, the list z will have the list of generators of 
the successive residue fields to be built, and the list Y will be the list of indeter- 
minates generating the corresponding polynomial rings. At this moment t is the 
order one type given by the lift 0i — x + 1 S Z[x] of y + l. We proceed to enlarge it, 
constructing a polynomial (/)2 £ Z[a;] whose 0i-Newton polygon has a unique side 
with slope 1/2 and residual polynomial j/^ + y + 1: 

> EnlargeType(l,2,y~2+y+l, t, Y, z) ; 

The polynomial 02 is 

> t[2]~Phi; 

x~4 + 4*x'3 + 8*x~2 + 8*x + 7 

We go one step further, and enlarge the type with a polynomial (\)-i £ Z[x] 
having a one-sided 02-Newton polygon, with slope 2/3 and residual polynomial 

y[2]2 + z[2]y[2] + z[l], where Ft = F2(z[2]) = F2[y]/(y2 + y + 1), and z\\\ = 1. 

> EnlargeType(2,3,Y[2]'2+z[2]*Y[2]+z[l] , t, Y, z) ; 

> t [3] ~ Phi ; 

x~24 + 24*x'23 + 288*x'22 + 2288*x~21 + 13482*x'20 + 
62664*x-19 + 238736*x-18 + 765072*x-17 + 2100447*x-16 + 
5005808*x-15 + 10455232*x-14 + 19267808*x-13 + 31474060*x-12 
+ 45694864*x-ll + 59023136*x-10 + 67789216*x-9 + 69069935*x-8 
+ 62166904*x-7 + 49095584*x-6 + 33671984*x-5 + 19743658*x-4 + 
9650728*x-3 + 3776208*x-2 + 1092816*x + 208081 

Let us build a second type: 

> tt,YY,zz:=InitializeType(2,y~2+y+l) ; 

> EnlargeType(3,2,y+l,--tt,--YY,--zz) ; 

> EnlargeType (3 , 1 , YY [2] ~2+YY [2] +zz [1] , ~tt , ~YY , --zz) ; 

> EnlargeType (1, 3, YY [2] +ZZ [2] ,--tt ,--YY,~zz) ; 

> tt[4] ~Phi; 

x~24 + 12*x"23 + 102*x'22 + 616*x~21 + 3045*x~20 + . . . 

We can combine any number of types to produce a polynomial whose family of 
attached types will be the given ones. 

> pol :=CombineTypes( [t ,tt] ) ; 

> pol; 

x~48 + 36*x"47 + 678*x"46 + 8808*x~45 +. . . 

You can check that the polynomial pol has the two types t,tt at the prime 
p = 2, and follow the flux of the algorithm using the flag montestalk to set the 
level of verbosity of the package: 
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> SetVerbose("moiitestalk",4) ; 

> K<w>:=NumberField(pol) ; 

> Factorizatioii(ideal(K,2)) ; 

Analyzing irreducible factor modulo p: YO + 1 



Analyzing type of order 1 

Phir= X + 1 

Sides of Newton polygon: [ 

[ -1/2, 0, 12, 24, ] 

] 



Analyzing side [ -1/2, 0, 12, 24, ] 
Slope: -1/2 
Origin: ( , 12 ) 
End point: ( 24 , ) 

Monic Residual Polynomial= Y0~12 + 

Y0~10 + YO'6 + Y0~2 + 1 
Factors of R.P.:= [ 
<Y0~2 + YO + 1, 6> 
] 

In some sense, these constructions lead to number fields specially ill posed for 
Montes algorithm, since the order of the types dominates the complexity of the 
computations. The relation: 

ttt 
deg(ir)= ^ l[eifl 

teT(/,p) fc=o 

between the degree of a number field and certain entries e and / (partial ramification 
indices and residual degrees) of the types attached to the prime ideals lying over the 
prime number p, shows that the growth of the order jjt of these types is logarithmic 
with respect to the degree. 

3. Description of the package 

3.1. New data types and attributes for ideals and number fields. The cen- 
tre of the attention of the package are number fields and its fractional ideals. Since 
computations tend to be quite time-consuming, the important information obtained 
from any operation is saved. This implies the introduction of new attributes for 
number fields, and new data types for handling their fractional ideals. 
The new attributes for number fields are: 

• FactorizedPrimes: List containing the rational primes already decom- 
posed in K. 

• Primeldeals: An associative array, with one entry per each rational prime 
p in the FactorizedPrimes; it contains the prime ideals in Zk dividing p. 

• Locallndex: Associative array which contains, for every rational prime 
in FactorizedPrimes, the local exponent of the index of the polynomial 
defining the number field. 
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• Discriminant The discriminant of the number field. 

• FactorizedDiscriminant: The factorization of the discriminant of the 
number field. 

• pIntegralBasis: An associative array, which for every rational prime p 
in FactorizedPrimes stores a p-integral basis of the maximal order of the 
field. 

• IntegralBasis: A list with a basis of the maximal order of the fields. 

These attributes are assigned and expanded when the corresponding objects are 
computed. For instance, an integral basis of K is not computed unless the function 
IntegralBasis is used, but then the last four attributes are set, and they are never 
recomputed. 

The fractional ideals are handled by means of two different record structures: the 
PrimeldealRecord and the IdealRecord. The first one is designed to contain the 
essential information produced by Montes algorithm for every prime ideal, while 
the second one is used for general ideals. Their concrete fields are: 

PrimeldealRecord : =recf ormat< 

Parent: // Number field where the prime ideal lies. 

Pol: // Polynomial defining the number field. 

Prime: // Rational prime p below the prime ideal. 

e: // Ramification index of the prime ideal. 

f : // Residual degree of the prime ideal. 

Generator: // The ideal is generated by Prime and this element. 

Position: Position of the prime ideal in the list 

K'Primeldeals [p] . 
Type: // Type corresponding to the prime ideal. 

>; 

IdealRecord : =recf ormat< 

Parent: Number field where the ideal lies. 

Generators: // Generators given in the definition of the ideal. 

Norm: // Norm of the ideal. 

Radical: // List of rational primes dividing the norm 

of the ideal. 
Islntegral: // Flag to describe the integrality of the ideal. 
IsPrime: // Flag to describe the primality of the ideal. 
Factorization: // List with the prime ideals dividing the 

ideal, and the corresponding exponents. 
FactorizationString: // Factorization of the ideal in 

pretty format . 
IntegerGenerator : // The ideal is generated by it and the 

element in the next field. 
Generator: // The ideal is generated by it and the element 

in the previous field. 

>; 

The fields in these structures are assigned only when its computation is trivial or 
strictly necessary. Once ideals or prime ideals are created, they can be manipulated 
using the standard notation of Magma for ideal arithmetic. 
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The trivial ideals (0) and Z/^ : an ideal is the zero ideal if its Generators field 
is assigned, and all its elements are zero; and an ideal is the full ring of integers if 
its Factorization field is assigned and is an empty list. One must check whether 
these fields are assigned, since the ideal may be the result of an operation between 
other ideals, which does not assign any of these fields. The functions IsZero tests 
whether an ideal is zero. 



3.2. Main functions. The essential tasks of the package is the creation and fac- 
torization of fractional ideals in number fields; they are performed by means of the 
functions ideal and Factorization. 

The purpose of the function ideal is the creation of a general fractional ideal, 
given a generator or a list of generators. Accordingly, it admits the two forms: 

ideal (K, alphal ); 

ideal (K , [alphal , . . . , alphar] ) ; 

where K is a number field and alphal , . . . , alphar are elements of K. The function 
returns an IdealRecord structure, with the fields Parent and Generators assigned. 
Once we have created a fractional ideal I, we can call the extended version of 
Factorization to find its decomposition into prime ideals. We can pass the ideal 
as a direct argument to the function or as a reference argument: 

Factorization(l) ; 
Factorization(~l) ; 

They do the same computations, but the procedural form assigns the fields 
Factorization, FactorizationString and IsPrime in the record given by 1 (if 
they are not yet assigned) and returns nothing, while the first form returns the value 
of the field Factorization. The behaviour of the functions depend on the data 
known for the ideal: If the field Factorization is assigned, the function returns 
its value. Otherwise, the greatest common divisors of the norms of the generators 
of the ideal is computed and factored. Every prime number dividing this gcd is 
decomposed in the number field (if it is not already decomposed) and the valuation 
of the given ideal at the prime ideals above these prime numbers are computed with 
the PValuation function. 

The decomposition of a rational prime in a number field is performed by the 
Montes function; the command 

Montes(p,K) ; 

applies Montes algorithm to the polynomial defining the number field K and the 
rational prime p, to find the prime ideal decomposition of p in K. It appends p to the 
list K~ Factor izedPrimes, and fills the entries corresponding to p in the associative 
arrays K~Locallndex, and K~Primeldeals. 

A key tool for the factorization is the PValuation function. It accepts an element 
of a number field and a prime ideal in the field (given as a PrimeldealRecord or 
an IdealRecord) and returns the valuation of the ideal generated by the element 
at the given prime. 

Since the discriminant of a number field is a very important invariant, the package 
includes the function TrueDiscriminant to compute it. Note that it has to factor 
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the discriminant of the polynomial defining the number field, which may be quite 
time-consuming or even impossible. 

3.3. Ideal arithmetic. We can add (+), multiply (*) or divide (/) two ideals and 
raise a non-zero ideal to any integral power (') using the standard operators. The 
mantra for the design of the arithmetic of ideals is "don't waste time in unnecessary 
hard computations". Thus, for instance, the addition of two ideals simply joins 
its list of generators and only if the factorization of both ideals is known also 
finds the factorization of the result. The multiplication of two ideals combines 
the factorizations of the ideals (computing them if necessary), but does not find 
generators for the product. 

The ideals can also be compared in the usual way: the operators eq and subset 
have extended versions for such task. They simply compare the factorizations of 
the ideals to be compared. The norm of an ideal and the rational primes dividing 
it are given by the functions Norm, and RationalRadical. 

3.4. Generators. It is well known that any fractional ideal in a number field can 
be generated by a rational integer and a second integral element in the field. The 
usual Magma function TwoElement has been extended to apply to our new ideals. 
It admits two forms 

TwoElement (I) ; 
TwoElement (~I) ; 

They perform the same computation, but the procedural version also assigns 
the fields Generators, IntegralGenerator and Generator of the given ideal. The 
factorization of the ideal is computed if it is not known. 

3.5. Congruences. Chinese remainder theorem. The package provides fa- 
cilities for congruential computations for prime ideals. Essentially, the package 
provides the reduction map Zx — > '^k/P from the maximal order of the number 
field K to the residue field of any prime ideal p. The residue field of a prime ideal 
is accessed through the command: 

ResidueField(P) ; 

The reduction mod p of a p-integral element a € K is computed with 

alpha mod P; 

An element z e Fp in the residue field of p may be lifted to an integral element 
in Zif with the command 

Lift(z,P) ; 

The Magma ChineseRemainderTheorem function is extended in the package to 
work with the new ideal records. Its syntax is 

CRT([alphal, . . . ,alphar] , [idl, . . . ,idr]) ; 

where alphal , . . . , alphar are integral elements of a number field and idl , . . . , 
idr are integral ideals (given as IdealRecord's or as PrimeldealRecord's). The 
function returns an element alpha in the same number field, such that alpha-alphak 
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is divisible by idk. Note that the function does not test the integraUty of the given 
elements. 

3.6. Integral Basis. A conjectural algorithm for the computation of integral bases 
of number fields can be derived from Montes algorithm. It has been developed in 
[GMNOQa] . The algorithm always returns a correct answer, since it checks its 
validity; in case of fail, the user is warned. Extensive tests of the algorithm have 
been done, without no counterexample found. 

The functions included in the package for this task are: 

pIntegralBasis(K,p) ; 
SIntegralBasis(K, [pi , . . . ,pr] ) ; 
IntegralBasis(K) ; 

The basic one is pIntegralBasis, which takes a number field K and a prime 
number p and computes a p-integral basis of the maximal order Zk; the basis is 
saved in the corresponding entry of the associative array K~ pIntegralBasis. The 
prime p is decomposed in K only if this has not been done yet. 

The SlntegralBasis functions takes a list of primes as second input, and returns 
a local integral basis for this set of primes. It calls pIntegralBasis for every given 
prime, and glues the local bases with a Chinese Remainder computation. 

Finally, the function IntegralBasis computes a Z-basis of the maximal order 
of the given number field. It simply calls SlntegralBasis with the list of primes 
dividing the discriminant of the polynomial defining the number field. 

3.7. Construction of hard number fields. In order to check the main functions 
of the package we have developed a number of functions to create examples of 
number fields with prescribed behaviour of the Montes algorithm. They allow to 
create and combines types from its essential data. The basic functions for this task 
are: 

InitializeType(p,psi) ; 
EnlargeTypeC h,e,psi, t, Y, z) ; 

A detailed example of its use has been presented in section 12.51 Remind that 
we can use the function CombineTypes to combine different types with the same 
rational prime. The function 

CombinePolynomialsWithDif f erentPrimesCf l,pl ,f 2,p2,k) 

takes two polynomials f 1 , f 2 of the same degree and two prime numbers pi , p2, 
and finds a polynomial / of the same common degree such that f = f (mod pi*^). 
It is useful for the construction of examples of number fields with prescribed types 
at different primes. 

Random examples can be easily produced using the functions: 

CreateRandomType(p,r) ; 

Great eRandomMultipleTypePolynomial(p,k,r, s) ; 

The first one generates a random type of order r for the prime number p, while 
the second one computes a random polynomial which has k different types of order 
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r for the prime p. This polynomial is obtained as the product of all the last (jj- 
polynomials of the randomly created types, plus an integral multiple of p~s to 
ensure its irreducibility. 

Since the examples produced by these functions use to be enormous, some func- 
tions to write them in a format proper for Magma or TeX are included. They 
are ExpandTeX, ExpandPhiTeX, ExpandAllPhiTeX and ExpandMagma. Their self- 
explanatory name makes unnecessary any further comment. 
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